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Abstract 

We study the production of hadrons in Au-|-Au collisions in the region 0.6 < xp < 

1.2, which we refer to as the trans-fragmentation region (TFR), since it corresponds 
roughly to r/' > 0, where rj' = r/ —T/beamj depending on pp- We show how hadrons can be 
produced in that region when the hadronization process is parton recombination. The 
inclusive x-distributions for proton and pion production are calculated with momentum 
degradation taken into acconnt. The results show that the proton yield is significantly 
higher than that of the pions in the TFR. Without particle identification the existing 
data cannot be used for comparison with our result on the p/vr ratio. Withont pp 
determination it is not feasible to relate the x distribution to the experimental r]' 
distribution. Nevertheless, on theoretical grounds we have shown why the production 
of hadrons in the TFR is not forbidden by momentum conservation. 

PACS numbers: 25.75.Dw 
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1 Introduction 


One of the interesting features of hadron production in heavy-ion collisions (HIC) is the 
discovery of the scaling property of the pseudorapidity distribution in the fragmentation 
region over an order of magnitude of variation in the collision energy [1, 2], Such a property 
was hypothesized in [3] for pp collision, and was referred to as limiting fragmentation. The 
independence on energy, however, does not imply independence on the collision system. The 
pseudorapidity distributions in the fragmentation region show dehnitive dependence on the 
centrality of Au-|-Au collisions for any hxed energy [2, 4]. Indeed, one cannot expect on 
theoretical grounds that the properties of hadron production in that region should be the 
same in pp and AA collisions. Whereas in pp collisions no particle can be produced with a 
rapidity greater than the beam rapidity, that is not the case in A A collisions. In fact, in terms 
of the variable p' = rj — yheam, where p is the pseudorapidity and |/beam the beam rapidity, 
PHOBOS data indicate that the charge-particle distribution dNch_/dp' does not vanish in the 
p' > 0 region [2, 5]. We shall refer to that region as the trans-fragmentation region (TFR). 
In this paper we describe the physics of hadron production in the TFR in the framework of 
the recombination model. 

Limiting fragmentation is a natural consequence of any formalism that uses momentum 
fractions of partons as the essential variables to describe hadronization since c.m. energy 
does not appear explicitly. However, particle production in the TFR is intriguing because 
it seems to violate momentum conservation, as it certainly would in pp collision. In AA 
collisions there are complications due to the fragments of the non-interacting spectators that 
can get into the detectors at small angles. For peripheral collisions at relatively low energy 
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= 19.6 GeV) dNch/dr]' seems to approach a constant value as 77 ' —> 2 [2, 5]. Such effects 
seem to diminish at higher ^/s. Let us put aside such issues, since our interest here is in the 
hadronization of the interacting part of nuclear collisions. For that, there is at present no 
data in the TFR for y/s = 200 GeV, say. Nevertheless, it is an interesting and important 
question to ask whether there exists theoretical reasons for hadrons to be produced in the 
TFR. 

Since the problem deals with \ow-pT physics, one cannot make use of pQGD with any 
degree of conhdence. Nevertheless, if one takes the point of view that hadrons are produced in 
the fragmentations region (FR) by parton fragmentation, such as in the dual parton model 
[ 6 ], then no hadrons can appear in the TFR, since all partons have momentum fractions 
a; < 1. On the other hand, in parton recombination the momentum fractions add and can 
result in a hadron momentum fraction greater than 1 , provided that the constituents come 
from different nucleons in the colliding system, a condition that is readily satisfied in AA 
collisions. It is with that possibility in mind that we study in detail the problem of hadron 
production in the TFR in the recombination model (RM) [7, 8 , 9]. 

Nuclear collisions have the complication of momentum degradation of constituents travers¬ 
ing nuclear matter. We have investigated the degradation effect in pA collisions, treating the 
constituents in terms of valons [ 8 , 10 ], and found good agreement with the data on “baryon 
stopping” in the FR [11, 12]. Here, in AA collisions the medium is dense, and our focus is 
in the TFR. Nevertheless, similar formalism will be used to take the degradation effect into 
account. 

Since hadronization is a problem that involves the momenta of the consituents in an 
essential way (as does momentum degradation), we shall be working with the momentum 
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fraction variable x, instead of the pseudorapidity variable rj. In terms of x the TFR is more 
precisely a; > 1, and the FR is for a; < 1, but above, say, 0.2. These regions do not map 
isomorphically to regions in 77 , since pt is involved in the dehnition of angle 6 . Although 
?7 is a more convenient variable for experimental detection, we shall work entirely in the x 
variables. The mismatch between theory and experiment in that respect can be overcome 
only when the pt values of the detected particles are determined. Until then, we cannot 
compare our predictions with any data. In our treatment the physics of hadronization on 
both sides of a; = 1 is continuous, so we shall calculate the spectra in the region 0.6 < a; < 1.2 
which we broadly refer to as TFR, as it roughly corresponds to the 77' > 0 region (although 
in principle, there is no upper limit on x). 

Apart from details our main qualitative prediction is that protons dominate the charge 
particles detected in the TFR. Such a prediction should be easier to conhrm or falsify than 
the verihcation of our results on the inclusive distributions in x. It is hoped that this paper 
will stimulate an effort to improve particle idenhcation in the TFR. 

2 Preliminary Considerations 

Let us begin with the kinematics x relevant for the TFR in HIC at high energy. For hadrons 
detected at small angle 6 relative to the beam axis, we may approximate tan 0/2 by pt/^pl, 
where pt and pl are the transverse and longitudinal momenta, respectively, of the produced 
hadron. Thus the pseudorappidity is 77 = \n{ 2 pL/pT), while the beam rapidity is 7 /beam = 
\n{^/s/nip), where nip is the proton mass. With 77 ' being the shifted pseudorapidity, 77 ' = 
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V — l/beam, we have for Feynman x, defined by a; = 2pl/ y/s, 



The mapping between rj' and x therefore depends on px- If in the forward region {px) < rup, 
then the p' > 0 region corresponds to a range in x that straddles a; = 1. For that reason we 
study the hadron spectra in the range 0.6 < a; < 1.2 as a representative of the TFR. 

Next we consider the geometrical aspect of nuclear collisions. For convenience in iden¬ 
tifying the forward direction, let us consider AB collision, and we are interested in the 
production of hadrons in the TFR of A. In the Glauber model we write for the cross section 
of AB collision in the form 


= I d^g^^ib), 


( 2 ) 


where 


+ J d^sTeis) [l - . (3) 

Ta{s) is the thickness function normalized to A, i.e., 

Ta{s) = ^ j dzp{s,z), J d^sTA{s) = A, (4) 

p being the nuclear density normalized to 1; cr in Eq. (3) is the inelastic nucleon-nucleon 
cross section. 

To calculate the average number of wounded nucleons in R, we consider only the first 
term in Eq. (3) and call it g^^^^{b). Dehne to be the probability of A having u 
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collisions in B at impact parameter b, so that 




Then we have 




where is the corresponding probability in pB collision [13] 


K^{b) = [(TTB{b)Y exp [-aTB{b)] , 


from which one can recover the necessary condition 


^7rf(6) = l-exp [-aTB 


as required by Eqs. (5) and (6). Moreover, we obtain 


^ = aTB{b) 


Thus for pB collision the average number of wounded nucleons in i? at 6 is 


_ aTBjb) 


pP^{b) = 


1-exp[-aTB(6)] ‘ 


Now, returning to AB collision the average number of wounded nucleons in B is 


_ El, aSdh T4s)TB(\g- 6|) 




gWB(b) 


Since the number of binary collision, N^^{b), is the numerator of the last expression in Eq. 
(11), and the number of participants, is g^^{b), we have for A A collisions 

N^^A{b)/2 ■ 
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For Au+Au collisions at 200 GeV, the tabulated values of A^^coii and Apart are 1065 and 351 
(at 0-5% centrality), and are 220 and 114 (at 30-40%) [14], Hence we get 

pAM = 6,1 (0 - 5%) , 

3.9 (30 - 40%) . (13) 

The impact parameters that correspond to the two centrality bins can be calculated from 
the overlap function 

TAAih) = I d^sTA{s)TA{\s-b\) . (14) 

Using the simplihed form for p{s, z) with uniform density in Eq. (4) results in a distribution 
for TAA{b) that is slightly higher than the tabulated values for various centrality bins given 
in [14], Nevertheless, from the shape of the distribtion the corresponding values of b can 
reasonably be set at 

6 = Ifm (0-5%) 

8 fm (30 - 40%) . (15) 

The consideration of nuclear geometry and the associated wounded nucleons will become 
important in the following when the momentum degradation effect is to be taken into account. 

3 Hadron Production at Large pi in the Recombina¬ 
tion Model 

In pp collisions no particle can be produced with pi > -/s/2. In AB nuclear collisions one 
would initially expect the same to be true also. However, n nucleons in A have a combined 
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momentum of n^/2, which can make possible a particle produced with > V^/2 without 
violating momentum conservation, if a coherence effect is at work. That would be the case 
if the valence quarks in three nucleons in A, each with momentum fraction Xi > 1/3, say, 
recombine to form a nucleon, whose Feynman x can then exceed 1. It is also possible for a 
pion to be produced with a; > 1, but since an antiquark in the sea (or a gluon converted to 
qq) is needed with large enough Xi to recombine with a valence quark of another nucleon, 
the probability is much lower. This type of consideration need not be restricted to the a; > 1 
region. Even for x > 0.6, which we broadly refer to as TFR, the recombination of partons 
from different nucleons in A will dominate over those processes where the partons are from 
the same nucleon, as in pB collisions. That dominance is over and above the shifted peak 
in rapidity due to “baryon stopping” simply on the basis of extra momentum availability in 
AB collisions. Thus even without detail calculations we can predict that p/tt ratio is large 
in the TFR in AB collision. Similarly, A/K ratio is also large. 

Hadron production at low pt and large pl in hadronic collisions has been treated in the 
RM in good agreement with data [7, 8, 9]. The extension now to nuclear {AB) collisions 
with emphasis in the TFR region has the same basic recombination formula 



= X 


dx 


dxi dx2 dx^ TpAB/ \ D r 

- Fuud{Xl,X2, X3)Rp{Xi, X2, Xs, X 

Xi X2 Xs 


(16) 


for the production of proton. Rp is the recombination function (RF) that has been studied 
in the framework of the valon model for the nucleon structure [8, 10], by relating Rp to the 
valon distribution Gp in the proton [15] 




where 


Gp{yi, 1/2, i/s) = 9 p{yiy 2 Ty^^{yi +1/2 +1/3 - 1), (I8) 

Qp = \B{cx + 1, Q; + /3 + 2 ) 5 ( 0 ; + 1, /3 + 1)] , (19) 

and Qst is the statistical factor 1/6 [15]. It is by successfully htting the CTEQ parton 
distribution functions at low that the parameters a and f3 are determined to be [10] 

a = 1.75, f3 = 1.05. (20) 

The 6 function in Eq. (18) enforces the momentum sum for recombination: = x. 

The key quantity in Eq. (16) is the 3-quark distribution F^J^{xi,X 2 ,X 3 ). If it were 
like in pB collision, then the quarks would all originate from the projectile, and the Xi in 
^ 2 , X 3 ) would satisfy J2iXi < 1. However, in the TER of AB collisions we consider the 
dominant component where each quark is from a separator nucleon in the same longitudinal 
tube at distance s from the center of A, so that we can write in the factorizable form 

= F-(x,)F-(x2)F^(xs) ( 21 ) 

with the labels s and b suppressed. On the LHS of Eq. (21) we use (3) instead of A in 
the superscript to emphasize that only 3 nucleons in A are considered. How to generalize 
from 3 to H will be discussed in Sec. IV. The effect of momentum degradation due to the 
passage through B will be considered in the next section. The crucial point here is that the 
variables Xi in Eq. (21) are independent of one another, so the integrals in Eq. (16) are from 
0 to 1 for each Xi. Thus the maximum possible a; is 3, well beyond the conventional FR. 
It is this unconventional possibility of producing a proton in the TER that motivates our 
investigation here. 
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In Eq. (21) we have used the superscripts u and d to denote the flavors of quarks that 
are to recombine in Eq. (16) where the RE is given by Eqs. (17) and (18) with yi and y 2 
referring to the u quark and y^ to the d quark. However, if the projectile A is an isoscalar 
which we shall assume, then at every impact parameter s there are equal numbers of protons 
and neutrons, so = Fd. In the valon model for pB collisions the quark distributions are 
111 ] 

F^'‘(x,) = £dy'G',(y')K(A , (22) 

where the valon distribution G'^{y') differs from G{y) due to momentum degradation, as will 
be discussed in detail in the following section. K{z) is the quark distribution in a valon. 
Since our valon distribution is flavor independent, K{z) consists of both valence and sea 
quarks [10] 

K{z) = K^sG) + L{z) , (23) 

where 

Kns{z) = z^^G-zf/B{a,h + l) , a = 0.35, 6 =-0.61 . (24) 

These parameters a and b are determined from the moments of K^s given in [10]. For the 
sea quark distribution L{z) we have to do two things: first, we have to average the favored 
(Lf) and unfavored (L^) distributions given [10]; second, we readjust the normalization to 
saturate the sea. Lf{z) and Lu{z) have been determined for the proton sea to distinguish, 
for example, a u quark in a U valon (favored) from a d quark in a U valon (or a u quark in 
a D valon, unfavored). They are, for Q = 1 GeV/c, [10] 

in Lf{z) = -2M + 0.08t-10AG-6G 
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In Lu{z) 


-2.92 + 4.0t-5.95t2-1.4t3 , 


(25) 


where t = — \n{l — z). We define 

L,(z) = i [Lf{z) + L^(z)] . (26) 

This distribution does not include the conversion of gluons to the sea quarks, a process that 
we must consider in order to account for the hadronization of all partons, including gluons 
[8, 9]. Thus to saturate the sea, we renormalize as follows 

L',(z) = ZL,(z) , (27) 

where Z is determined by solving the two algebraic equations that express the momentum 
conservation in the u, d, s sections in terms of the second moments [9] 

K^si2) + 2 [2Z,(2) + L,(2)] + L,(2) = 1 (28) 

i^A7s(2) + 2[2L;(2)+L((2)] =1 (29) 

Z = 1 + r - ^ . (30) 

2[2L,(2)+L,(2) 

In Eq. (29) we have assumed that the s-quark sea is enhanced also in nuclear collisions, 
unlike the case of hadronic collision where the gluons are convected to the light-quark sector 
only through g ^ qq [9]. In getting Eq. (30) we have set = ZLg also. Since the second 
moments in Eq. (30) are tabulated in [9], we obtain 

Z = 3.42 . (31) 

In the following we shall consider L'^{z) only whenever the light quarks in the sea are needed, 
including the L{z) term in Eq. (23). 
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For pion production the physical content of the calculational procedure is the same as 
for proton production, except that it is a qq recombination. Thus, as in Eq. (16), (21) and 
(22) we have 

H^^{x) = = f ——F^^^{xi,X 2 )Rn{xi,X 2 , x) , (32) 

ax J Xi X 2 

= , (33) 

F|(i2)= /‘%'G„(!,')q(^) . (34) 

Jx 2 \y / 

where Fp is obtained from the saturated sea. The recombination function for pion is [8, 9] 

x) = ^6(^ + ^-i) (35) 

X^ \ X X ) 

In an isosymmetric collision system we need not distinguish the charge states. A and K 
production can similarly be considered. 

4 Momentum Degradation 

In the preceding section we have described the quark distribution in Eq. (22) as a convolution 
of the valon distribution and the quark distribution in a valon, but we have not specihed the 
former. In a pp collision the valon distribution G{y), y beginning the momentum fraction, 
has been identihed as that of the free proton on the basis that at low pt and large pi the fast 
partons in the forward direction are unaffected by the opposite-going partons due to the lack 
of long-range correlation in rapidity. The valon model connects the bound-state problem of 
a static proton (in terms of constituent quarks) with the structure problem of a proton in 
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collision (in terms of partons) [8]. In pA collision the effect of momentum degradation in 
the passage of the projectile through the nuclear target is applied to the valons in Ref. [11], 
where it is shown not only how baryon stopping can be obtained in agreement with data, 
but also how pion production in the FR can be determined. We now extend the treatment 
to the TFR in the AB collisions. 

In a free proton the single-valon U and D distributions are obtained from the exclusive 
distribution given in Eq. (18) by integration 

jj r^-yi ri-yi-y2 

G {yi) = dy 2 c^l/s Gp(|/i, 1 / 2 , 1 / 3 ) 

JO JO 

= (7pR(a + l,/3+l)|/“(l-|/i)“+^+i (36) 

„ r^-ys r^-yi-ys 

G (ys) = dyi d |/2 Gp(|/i, 1 / 2 , 1 / 3 ) 

JO JO 

= gpBia + 1, a + l)yi{l - 1 / 3 )'“+' • (37) 

For isoscalar nuclei we take the average 

G{y) = 1 [G‘'(!/) + G°(!/)] (38) 

before the nucleon traverses the nuclear medium. We note that G{y) is not an invariant 
distribution, but yG{y) is. It is normalized by 

^ dyG{y) = J dyidy 2 dy 3 G{yi,y 2 ,y 3 ) = 1 , (39) 

i.e., the probability for the proton to consist of three and only three valons is 1. 

Suppose now that a nucleon in R at a fixed impact parameter s collides with 9 wounded 
nucleons in B on the average. Since |s — b\ can be almost as large as the radius R, it is 
possible for 9{b, s) to be very small. That is the geometrical situation where several nucleons 
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in A can contribute to the TFR without much momentum degradation. When we consider 
fluctuations of u from 9, we must include the possibility of u being 0. Thus we use the 
Poisson distribution 

P,{u) = (40) 

with normalization defined by summation from z/ = 0 

OO 

= (41) 

i /=0 

If a valon loses a momentum fraction 1 — k at each collision, then after u collisions the 
modified valon distribution is 

v'GUy') = l^^dyG{y)s(^-A , (42) 

from which follows 

GUy') = . (43) 

The parameter k for the reduced momentum fraction is unknown, since it should not be 
inferred from pA collision, which is for a cold nuclear target. In AB collision a tube in 
A contains many nucleons which cannot all be treated as if each of them collides with u 
nucleons in a cold nucleus B. Indeed, it is hard to assess the state of B when the back part 
of the tube traverses the medium. We can at best use an adjustable parameter k to describe 
in some average sense what contributes to the TFR. From Eq. (43) we get after 9 collisions 
on the average 

OO 

G'M) = i: Gliy’)P,W) . (44) 

i /=0 
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To identify v with the quantity expressed by Eq. (10) is based on the assumption that each 
valon experiences the same average number of collisions as the parent nucleon does. 

Gp (?/'), as given in Eq. (44) is the modihed valon distribution that should be used in Eqs. 
(22) and (34). That takes care of the nuclear effect. What remains is the calculation of the 
p and TT distributions in the TER. 

5 Proton and Pion Distributions 

To calculate the hadron distributions in x, we return to Eq. (16) as a general formula, which 
needs, however, some more elaboration to account for the size of A. In Eq. (21) we show 

/'Q\ O 

the factorizable form of when there are only three nucleons in the projectile, each 

contributing a quark. Now, we specify the details of how to calculate that is called for 
in Eq. (16). 

In Sec. II we have in Eq. (7) the probability vrj)^(&) for a nucleon making u collisions in 
B at impact parameter b. Now, consider the same quantity in A and write 

^ WTAis)f exp [-aTAi-s)] (45) 

jj,'. 

for the probability that p nucleons in A colliding with a nucleon in B at impact parameter 
s in A. In place of Eq. (6) we now have for p nucleons in A colliding with v nucleons in B 
the probability 

n;!«(4) = /K|). (46) 

Since at least 3 nucleons in A are needed for our calculation of the proton distribution in 
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the TFR of A, as done in Sec. Ill, we must sum over fi starting with fi = 3 and get 


&) = / ^ £ 6, s) , 


where s) is what we have described in Eqs. (16), (21) and (22), i.e.. 




n 


u=i 


Rp{xi,X 2 ,X 3 ,x) , 


where z/ is given by Eq. (10) but at |s — 6| 


aTeds - b\) 


p = iyP^(\g-b\) = - 

1 - exp [-(tTb(|s - b\) 

The sum over jj, in Eq. (47) can be performed, yielding 


H^^(x.b) = 




( 47 ) 


(48) 


(49) 


(50) 


Note that with this formula we do not need the results on given in Eqs. (12) and 

(13). 

For pion production it is straightforward to modify Eq. (50) and get 


H^^(x,b) = 




(51) 


where 


Hl^^^{x,b.s) = [——F!j{xi)F^{x2)RTrixi,X2,x), (52) 

J Xi X2 

Since F^{x 2 ) is severely damped at large X 2 -, H:^^{x,b) is expected to be much more sup¬ 
pressed compared to Hp^{x,b) in the TFR. Nevertheless, H^^[x,b) need not vanish for 
a; > 1, unlike HP^{x,b). 

On the basis of Eqs. (50) and (51) we have calculated the proton and pion distributions 
for Au-|-Au collisions at 6 = 1 and 8 fm, corresponding to 0-5% and 30-40% centralities. 
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according to Eq. (15). We have used the approximation of uniform nuclear density with 
R = fm, and a = 41 mb. For the parameter k that represents the surviving valon 

momentum fraction after each collision, we have chosen two representative values, k = 0.8 
and 0.6, where k = 1 implies no momentum degradation. The results for the invariant 
distribution xdNp/dx for the proton, which is just Hp'^^'^{x,b), are shown in Fig. 1 for (a) 
& = 1 fm and (b) 6 = 8 fm. The solid (dashed) line corresponds to k = 0.8 (0.6). All four 
lines are nearly straight, i.e., exponential in x, smoothly throughout the TFR. There can 
be other contributions to the proton production in that region due to the recombination of 
quarks originating from one or two nucleons, but they are so small that we ignore them. 

Figure 1 shows that the x distributions are suppressed when there is more momentum 
degradation (smaller k), as is expected. In the case of pA collisions such a suppression would 
correspond to the qualitative notion of baryon stopping. But in AA collisions, instead of 
stopping, we have protons produced at a; > 1. Nevertheless, the overall normalization is 
lowered when there is more momentum degradation. Thus there are two features about 
the inclusive distributions of the proton: it extends smoothly beyond x = 1, and is more 
suppressed at lower k. The physical value of k, that corresponds to reality can be deter¬ 
mined only after data become available and put on plots such as Fig. 1. For the purpose of 
summarizing the behavior of the x distributions, we parametrize them in the form 

dN 

x—^ih) = exp [ho(&) — h\ih)x \, 0.6 < a; < 1.2 , (53) 

dx 

which hts the lines in Fig. 1 extremely well with the parameters /io(&) and hi{b) given in 
Table 1. 

It is of interest to know how the distribution depends on centrality at a given k. We plot 
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Table 1; Parameters /ig and h\ for proton 



0 

bo 

K = 0.6 


b = i 

00 

b = l 

00 

14 (b) 

10.58 

9.90 

7.03 

7.93 

im 

9.80 

8.18 

8.41 

8.04 


the distributions normalized by iVpart/2 in Fig. 2, which shows how they are suppressed as 
the collision changes from peripheral to central. At larger b the collisions on average have 
higher /i and lower u in Eq. (46), leading to more partons hadronizing in the TFR. 

In [2] there are data from Au+Au collisions at 200 GeV showing that the r]' distributions 
for 0-6% and 35-40% cross over at rj' ^ —2. According to Eq. (1) that would correspond to 
X ~ O.lSSpr/mp ~ 0.1, which is signihcantly outside the TFR. For —2 < 77 ' < 0 the data 
show the peripheral case higher than the central case, not unlike what we have in Fig. 2(a). 

For pion production Fig. 3 shows the rapid decline of the x distributions, as x is increased 
toward 1 and beyond. Since q distribution is suppressed at large X 2 , it is difficult for a pion 
to acquire enough momentum to go deep into the TFR. Nevertheless, the a; > 1 region is not 
forbidden. It is evident from Fig. 3 that there is no sensitive dependence on k. The reason is 
partly because there are only two G'p functions in Eqs. (33) and (34) for n, instead of three 
for p, but mostly because the RF for tt is broader than that for p. The valon distribution 
in TT is flat [9] (corresponding to pion being a tightly bound state), so the wider momentum 
spread allows the more degraded parton momenta to contribute to the formation of tt. In 
Fig. 4 we show the dependence on b. Compared to Fig. 2, the pions do not show as much 
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dependence as do the protons, although the vertical scales of the two figures are different 
and can lead to visual misreading. Quantitatively, we can £t the distributions by 


= exp ho(6) + h\{h)x — h 2 {h)x^ 


0.6 < a ; < 1.2 , 


(54) 


with the parameters given in Table II. Qualitatively, the pion distributions are orders of 


magnitude lower than the proton distributions. 


Table 2: Parameters hg, hi and for pion 



K = 0.8 

K = 0.6 


h = l 

00 

6 = 1 

00 

K 

1.96 

1.25 

0.0063 

-0.083 

K 

5.46 

7.39 

7.40 

8.32 

hi 

14.15 

14.79 

14.66 

14.87 


Apart from the details of the x distributions, our main prediction is that proton produc¬ 
tion dominates over pion production in the TFR. To give a visual impact of that dominance, 
we show in Figs. 5 and 6 the p/tt ratio. At a; ~ 1, the ratio is roughly 10^ for any combination 
of b and k. It is such a large ratio that particle identification in the TFR would be the most 
direct way to settle the question whether our hadronization scheme is in any way close to 
reality. 
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6 Conclusion 


We have investigated hadron production in the trans-fragmentation region. The overwhelm¬ 
ing feature of our result is that the p/vr ratio is extremely large, roughly 10^ at a; ~ 1. That 
feature is a direct consequence of parton recombination. For a proton to be produced at 
a; ~ 1, it is rather easy to hud three nucleons in A each contributing a quark at a;* ~ 1/3 
to form the proton. However, for a pion at a; ~ 1, it is hard to hnd any antiquark at, for 
example, a;* ~ 1/3 to help a quark at ~ 2/3 to make up the pion momentum. Quantitative 
value of the p/tt ratio can be obtained only after some rather involved calculations. But 
qualitatively to have a distribution in x that crosses the boundary at a; = 1 smoothly is 
possible only by parton recombination, since fragmentation would require all hadrons to be 
produced at a; < 1. 

Experimental data do show that particles can be produced at rj' > 0, at least at lower 
energies. At 200 GeV the data have stopped at p' ~ 0, but show no evidence of vanishing 
there [2]. It is unfortunate that we cannot compare our result in x to the data in p', since pt 
of the produced particles are unknown. Either particle identihcation or p^ determination, 
preferably both, would greatly help to relate theory and experiment. 

The importance of clarifying what happens in the TER is in the determination of whether 
there exists new physics in the FR and beyond. In our treatment of the problem we have 
considered only low-p^ physics, but extended to include recombination of quarks at medium 
Xi from different nucleons. If proven correct by data, it forms the basis from which to 
extend further to higher pt in the TFR. Then there should arise a competition between 
the enhancement effect studied here and the suppression effect found earlier in forward 
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production at intermediate pt in d+Au collisions [16]. The suppression can be due to either 
initial-state or hnal-state physics. Any future study of hadron production in a larger domain 
in Pt in the TFR will have to be consistent with the physics explored here at low pt- 
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Figure Captions 

Fig. 1. Proton distributions in the TFR for (a) & = 1 fm, and (b) 6 = 8 fm. 

Fig. 2. Same as in Fig. 1 but for (a) k, = 0.8, and (b) n = 0.6, where k, is the survival factor 
in momentum degradation. 

Fig. 3. Pion distributions in the TFR for (a) 6 = 1 fm, and (b) & = 8 fm. 

Fig. 4. Same as in Fig. 3 but for (a) k, = 0.8, and (b) k = 0.6. 

Fig. 5. p/t: ratio for (a) 6 = 1 fm, and (b) 6 = 8 fm. 

Fig. 6. Same as in Fig. 5 but for (a) k = 0.8, and (b) k = 0.6. 
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